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PREFACE 

This  Memorandum  was  prepared  for  the  Vela  Analysis  project  under 
the  Advanced  Research  Projects  Agency's  contract  with  RAND.  The  pro¬ 
ject  is  a  broad  study  of  systems  to  detect  above-ground  nuclear  bursts, 
and  of  necessity  includes  research  into  satellite  systems.  The  Memor¬ 
andum  should  assist  in  determining  whether  or  not  gravity  gradient 
attitude  stabilization  is  adequate  for  particular  missions  under 
consideration. 
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SUMMARY 

Satellite  attitude  stabilization  systems  utilizing  the  gradient  of 
the  earth's  gravitational  field  have  the  advantage  of  being  completely 
passive  (i.e.,  no  power  is  required).  However,  the  stabilizing  moments 
developed  are  extremely  small  and  as  a  result  such  systems  are  subject 
to  attitude  perturbations  from  sources  which  would  ordinarily  be  neglected 
in  the  evaluation  of  an  active  stabilization  system. 

This  Memorandum  investigates  certain  of  these  sources  of  attitude 
perturbation  including  micro-meteoroid  impact,  solar  radiation  pressure, 
station -keeping  propulsion,  orbital  eccentricity,  and  on-board  rotating 
machinery.  In  the  analysis,  formulas  are  developed  which  determine  the 
magnitude  of  the  pitch,  roll  and  yaw  perturbations  resulting  from  each 
of  the  above  sources.  Application  of  these  results  to  an  assumed  ve¬ 
hicle  configuration  shows  that  the  resulting  attitude  perturbations  can 
be  of  the  crder  of  several  degrees.  While  these  perturbations  do  not 
necessarily  rule  out  gravity  gradient  stabilization,  their  magnitudes 
for  a  given  vehicle  configuration  should  be  considered  for  their  com¬ 
patibility  with  the  mission  requirements. 
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I.  INTRODUCTION 

It  has  been  recognized  for  many  years  that  a  rigid  body  tends  to 
align  itself  with  its  axis  of  minimum  moment  of  inertia  in  the  direc¬ 
tion  of  the  gradient  of  a  gravitational  field.  Until  relatively  re¬ 
cent  times,  this  effect  was  primarily  of  interest  as  an  explanation 
of  the  fact  that  the  moon  maintains  the  same  face  toward  the  earth. 
However,  with  the  development  of  artificial  satellites,  it  was  sug¬ 
gested  (see  Ref.-  1)  that  this  same  effect  could  be  used  to  stabilize 
such  vehicles  so  that  the  orientation  relative  to  the  earth  remains 
the  same  at  all  times.  The  feasibilir  •  of  such  a  system  for  some 
applications  has  recently  been  demonstrated  (see  Ref.  2). 

The  use  of  the  gravity  gradient  to  achieve  attitude  stability  has 
the  attractive  feature  that  the  system  is  completely  passive  and  thus 
puts  limitation  on  vehicle  lifetime.  However,  since  the  stabilizing 
moments  developed  are  dependent  upon  the  nonuniformity  of  the  gravita¬ 
tional  field  over  the  dimensions  of  the  vehicle,  it  is  obvious  that 
these  moments  are  extremely  weak.  Thus,  a  relatively  small  disturbing 
moment  might  override  the  inherent  stability  of  the  vehicle. 

This  Memorandum  investigates  the  effect  of  several  such  sources  of 
perturbations  on  the  performance  of  a  gra/ity  gradient  stabilization 
system.  The  perturbations  considered  include  the  effects  of  auxiliary 
rotating  equipment,  orbital  eccentricity,  as  well  as  external  moments 
due  to  such  factors  as  radiation  pressure,  micro-meteoroid  impact,  and 
station-keeping  propulsion. 


II.  ANALYSIS 


STATEMENT  OF  HIE  PROBLEM 

For  the  purposes  of  this  Memorandum,  it  is  assumed  that  the  satellite 
vehicle  is  on  a  nearly  circular  orbit  and  its  attitude  is  specified  by 
the  coordinate  systems  shown  in  F^‘g.  1.  The  X,  Y,  Z  system  is  an  or¬ 
bital  coordinate  system  with  the  Y  axis  along  the  instantaneous  vertical 
and  the  Z  axis  along  the  normal  to  the  orbital  plane.  The  x,  y,  z 
system  is  a  sec  of  body  axes  fixed  in  the  vehicle  and  with  an  orientation 
relative  to  X,  Y,  Z  specified  by  the  angles  a,  B  and  i|f  in  pitch,  roll 
and  yaw  respectively.  By  a  suitable  choice  of  the  moments  of  inertia 
in  pitch,  roll  and  yaw,  the  vehicle  will  be  stable  when  a,  B  and  ijf  are 
all  zero.  The  problem  then  is  to  determine  how  much  this  stability  is 
disturbed  by  various  perturbing  effects. 

GENERAL  EQUATIONS  OF  MOTION 

The  general  equations  describing  the  motion  of  the  satellite  about 
its  center  of  mass  are  derived  in  Ref.  3  and  can  be  expressed  as 


3g  R2 
o  o 


b  c 

y  y 


0) 


(2) 


(3) 


where  uj^,  u> ,  0)^  are  the  components  of  the  vehicle  angular  velocity 
relative  to  inertial  space  and  are  given  by 


Fig.  1  — Coordinate  systems 
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a;  =-(&  +  &)  cos  3  sin  i'  4-  B  cos  A  (4) 

oi^  =  (<y  +  3)  sin  8  +  f  (5) 

tD  =  (cy  +  0)  cos  8  cos  \}f  4-  8  sin  \}?  (6) 

2 

while  a  ,  b  ,  c  are  the  direction  cosines  of  the  Y  axis  in  the  x,  y,  z 

y 1  y 7  y 

system  and  are  given  by 


a 
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sin  (y  cos  +  cos  jy  sin  8  sin  i]j 
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cos  a  cos  g 


=  sin  a  sin  A 

y  * 


cos  a  sin  8  cos  ^ 


(7) 

(8) 
(9) 


The  quantities  I  ,  I  and  I  are  the  principal  moments  of  inertia  about 
x  y  2 

the  roll,  yaw  and  pitch  axes,  while  and  are  the  corresponding 

external  moments  about  these  same  axes.  In  addition,  R  is  the  radius 

o 

of  the  earth,  gQ  is  the  gravitational  acceleration  at  the  earth's  sur¬ 
face,  r  is  the  radial  distance  from  the  earth's  center  to  the  center 
of  mass  if  the  vehicle  and  0  is  the  instantaneous  orbital  angular  rate. 
LINEARIZED  EQUATIONS  OF  MOTION 

Since  in  any  satisfactory  stabilization  system  large  angular  de¬ 
viations  in  attitude  cannot  be  tolerated,  it  is  permissible  to  linearize 
the  equations  of  motion  developed  above  by  the  following  substitutions 

u>  =  -  0f  +  £  (10) 

=  88  +  f  (11) 

<«z  =  9  +  o- 


(12) 
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while 
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=  a 


(13) 
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y 


=  1 


(14) 


c 

y 


8 


(15) 


In  addition,  for  a  nearly  circular  orbit,  the  orbital  angular  rate  is 


given  by 


o  g  R  o 
Zoo  2 
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=  tu 


(16) 


Substitution  of  Eqs .  (10)~(16)  in  Eqs.  (l)-(3)  gives 
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These  equations  can  be  further  simplified  by  defining  the  following 
parameters 

Rx  =  '\~L  (20) 
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R 
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R 
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(21) 


(22) 


Substitution  of  Eqs.  (20)-(22)  into  Eqs.  (17)-(18)  results  in 
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(25) 


It  should  be  noted  that  the  quantities  R  ,  R  and  R  are  not  inaepen- 

x  y  z 

dent  and  it  is  possible  to  express  R  in  the  form 

z 


R  = 


R  -  R 

x  y 

1  -  R  R 
x  y 


(26) 


In  addition,  there  are  certain  other  limitations  on  the  values  of  R 

x 

and  R  dictated  by  (1)  the  combinations  of  I  ,  I  and  I  which  can  be 
y  x  y  z 

realized  physically,  and  (2)  the  requirement  that  the  system  be  in¬ 
herently  stable  in  pitch,  yaw  and  roll.  These  limitations  and  the 
natural  frequencies  in  pitch,  yaw  and  roll  are  determined  in  the  Appen¬ 
dix  . 

PERTURBATION  EFFECTS 

In  this  section  the  effect  of  various  disturbances  on  a  gravity 
gradient  stabilized  satellite  vehicle  will  be  considered.  The  dis¬ 
turbances  included  are  (1)  an  arbitrary  impulse  moment  such  as  might 
result  from  micro -meteoroid  impact  or  from  station -keeping  propulsion, 
(2)  a  step  increase  in  the  external  moment  such  as  that  due  to  radiation 
pressure  or  a  steady -state  thrust,  (3)  eccentricity  of  the  orbit  and 
(4)  the  interaction  between  the  vel  icle  and  on-board  rotating  equip¬ 


ment. 
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Impulse  Moment 

Pitch  Excitation.  If  it  is  assumed  that  an  impulsive  moment  of 

magnitude  1^  is  applied  about  the  pitch  or  z  axis,  then  the  external 

moment  M  is  given  by 
z 

Mz  =  Ip6(0)  (27) 

where  6(0)  is  a  Dirac  delta  function  at  t  =  0. 

Equation  (25)  can  now  be  written  in  the  form 


•  •  2 

a  +  u)  a  =  —  6(0) 
a  z 


(28) 


from  which  the  pitch  response  is  obtained  as 


y  -  - -  sin  ti)  t 

I  Hi  Of 

z  y 


(29) 


Thus,  the  maximum  pitch  angle  is  given  by 
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(30) 


which  can  be  rewritten  by  means  of  Eqs .  (16)  and  (191)  in  the  form 


y  = 
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I  R 
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r  (1  -  R  R  ) 
c  x  y 

3g  (R  -  R  ) 
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(31) 


where  r  is  the  orbital  radius, 
c 

Roll  Excitation.  If  an  impulsive  moment  is  applied  about  the  roll 
or  x  axis  ,  the  value  of  M  is  given  by 

X 

M  =  I  6(0)  (32) 

x  R 

Since  the  roll  and  yaw  motions  are  coupled,  it  is  necessary  to 
solve  Eqs.  (23)  and  (24)  simultaneously.  These  equations  can  be 
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rewritten  in  the  form 
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where 
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?e  solution  to  Eqs .  (33)  and  (34)  is  given  by 
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1  COS  UD^t 
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(33) 

(34) 

(33) 

(36) 

(37) 

(38) 

(39) 

(40) 


where  uj^  and  are  t^e  two  characteristic  frequencies  of  the  yaw-roll 
response,  as  determined  in  the  Appendix. 

Since  both  the  roll  and  pitch  responses  are  made  up  of  the  sum  of 
two  sinusoidal  terms  of  frequencies,  uj^  and  the  maximum  possible 
value  would  be  the  sum  of  the  amplitudes  of  the  individual  terms . 
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(41) 


while  the  associated  yaw  response  has  a  maximum  amplitude  of 


21 
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It  should  be  noted  that  Eq .  (41)  is  valid  only  if  R  and  R  are 

x  y 

both  positive,  which  is  the  case  for  most  configurations  considered 
for  gravity  gradient  stabilisation.  If  R  and  R  were  negative,  the 
second  fraction  in  the  brackets  in  Eq.  (41)  would  have  a  negative 
sign. 

Yaw  Excitation.  If  an  impulsive  moment  is  applied  about  the  yaw 
or  y  axis,  the  resulting  value  of  is  given  by 


My  -  y(°; 


(43) 


and  Eqs .  (23)  and  (24)  have  the  form 


3  +  <Ug3  +  M'p'jf  =0 


(44) 


M-.B*  +  ijr  +  uA  =  6(0) 

XT 


(45) 


The  solution  for  the  resulting  yaw  and  roll  responses  is  given  by  the 
relations 
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And 
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V'h  -  »2> 


COS 
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(47) 


As  in  the  previous  case,  the  maximum  possible  amplitudes  in  yaw  and 
roll  are  given  by 
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As  in  the  case  tor  roll  excitation,  Eq .  (48)  applies  only  when  R  and 

X 

R  are  positive.  If  R  and  R  are  negative,  the  sign  of  the  second 

y  *  y 

fraction  in  brackets  reverses. 

It  should  be  noted  that  in  this  and  the  succeeding  analyses,  no 
damping  has  been  incorporated  in  the  control  system.  If  this  were  in¬ 
cluded,  the  transients  resulting  from  impulse  moments  would  eventually 
damp  out.  However,  the  expressions  given  for  the  maximum  values  of 
O',  8  and  4  are  a  good  measure  of  the  attitude  disturbances  which  might 
exist  during  the  early  part  of  the  transient  response  since  the  time 
ccistant  of  a  system  with  damping  is  of  the  order  of  several  orbital 
periods . 

Step  Moment 

Pitch  Excitation.  If  a  step  change  occurs  in  the  moment  applied 
about  the  pitch  axis  such  that 


M  -  0 
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(t  >  0) 


(50) 
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then  Eq .  (25)  can  be  written  as 
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tV  +  MJ  C*1 
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(51) 


which  has  a  solution  given  by 


a  - 


I  HI2 
z  ot 


(1 


cos  <J0  t) 
at 


(52) 


It  is  seen  that,  in  the  absence  of  damping,  the  pitch  angle  reaches 
a  maximum  of  ^ 
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If  damping  is  present,  then  the  steady-state  value  of  the  pitch  angle 
is  half  that  shown  above  or 
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Roll  Excitation.  If  a  step  change  occurs  in  the  moment  applied 
about  the  roll  axis,  Eqf .  (23)  and  (24)  representing  the  coupled  yaw- 
roll  response  can  be  written  as 
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The  resulting  yaw  and  roll  responses  are  given  by 
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An  examination  of  Eqs .  (57)  and  (58),  together  with  Eqs .  (21),  (36), 
(193)  and  (194),  shows  that  the  oscillatory  components  of  the  responses 
have  maximum  possible  amplitudes  given  by 
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In  addition  to  the  transient  terms  discussed  above,  the  roll  response, 
Eq .  (57),  also  has  a  steady-state  term  given  by  the  expression 
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It  is  seen  that  3  and  3  are  equal;  thus,  as  a  worst  c'.se  the 
m  ss  ’  ’ 

roll  angle  might  have  a  value  equal  to  twice  that  given  in  EH  .  (61) 
Yaw  Excitation.  If  a  step  change  in  the  yaw  moment  occurs, 

Eqs .  (23)  and  (24)  take  the  form 
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(63) 


By  analogy  with  the  previous  case,  the  yaw  and  roll  responses  can  be 
written  as 
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and 
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sin  aj  t 
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sin  u^t 


(65) 


The  maximum  possible  amplitudes  of  the  oscillatory  terms  in  Eqs.  (64) 
and  (6. i)  are  determined  as  in  the  previous  case,  and  are  of  the  form 
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In  this  case,  there  is  a  steady-state  deflection  in  yaw  given  by 
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Again,  the  values  of  $  and  j*  are  equal  and  a  yaw  angle  equal  to 


m 


ss 


their  sum  might  result. 

It  should  also  be  noted  that  if  the  pitch  and  roll  moments  of  inertia 
are  nearly  equal,  the  resulting  value  of  $  may  be  large. 

Orbital  Eccentricity 

In  the  development  of  the  linearized  equations  of  motion,  it  has 

been  assumed  that  the  orbital  angular  rate  0  is  a  constant,  to  .  How- 

o 

ever,  if  the  orbit  has  a  small  eccentricity,  @  is  no  longer  constant 
and  the  linearized  equations  obtained  by  the  substitution  of  Eqs. (10)-(15) 
into  Eqs.  (l)-(3)  have  the  form 


3  +  Rx  020  +  (Rx  -  1) 


-  „2 
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(69) 


|  0  i<  -  ( R  -  1 )  98  =  0 


(70) 


15 


3  g  R2 

0  +  a  =  -  - y~  \ot 


(71) 


The  orbital  motion  of  the  satellite  is  specified  by  the  relations 


i  8  R‘ 

1  o  o 


■55—  (1  +  e  cos  9) 


(72) 


re  =  c 


(73) 


where  C  is  a  measure  of  the  orbital  angular  momentum,  e  is  the  orbital 
o 

eccentricity  and  the  angle  0  is  measured  from  perigee. 

Pitch  Response.  By  means  of  Eq .  (73),  it  is  possible  to  change 
the  independent  variable  of  Eq .  (71)  from  t  to  0  so  that 
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=  0 


(74) 


Elimination  of  r  between  Eqs.  (72)  and  (74)  gives 

2 

+  3R  or  =  2e  sin  0  (75) 

d0 

where  second-order  terms  have  been  neglected. 

It  is  seen  that  the  eccentricity  appears  as  the  amplitude  of  a 
sinusoidal  driving  function  in  Eq .  (75)  and  the  resulting  pitch  re¬ 
sponse  is  given  by  the  relation 
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(76) 


The  maximum  possible  pitch  angle  is  expressed  as 
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If  the  system  has  damping,  the  second  tern  of  Eq .  (76)  will  eventually 
damp  out  and  the  resulting  pitch  amplitude  will  be  given  as 
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m  6.-1 
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2e 
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(78) 


An  examination  of  Eq .  (76)  shows  that  this  solution  is  not  valid  if 

is  equal  to  either  zero  or  one-third. 

The  first  of  these  cases  (R  =  0)  corresponds  to  a  configuration  which 

is  symmetrical  about  the  pitch  axis  and  thus  has  no  pitch  stability.  In 

the  limit  as  R  approaches  zero,  Eq .  (76)  reduces  to 
z 

Of  =  2e  |*0  -  sin  6]  (79) 


Since  this  solution  has  a  secular  term,  sooner  or  later  it  would 

violate  the  assumption  that  a  is  a  small  angle.  In  Fig.  3  of  the 

Appendix,  the  locus  for  this  case  is  the  contour  for  cr  equal  to  zero. 

O' 

In  the  second  case  (R  =  1/3)  the  term  on  the  right  in  Eq.  (75)  is 

z 

driving  the  system  at  resonance,  and  as  R  approaches  1/3,  Eq .  (76) 

z 

takes  the  form 

0/  =  e  [sin  9-0  cos  0]  (80) 

This  constitutes  a  divergent  oscillation  which  would,  in  time,  exceed 
the  small  angle  assumption  for  a.  This  case  corresponds  to  the  contour 


17 


for  m  /u>o  equal  to  unity'  in  Fig.  3  of  the  Appendix. 

Yaw-Roll  Response.  If  Eqs.  (b9)  aid  (70)  ore  expressed  with  8  as 
the  independent  variable  and  second-order  terms  are  neglected  as  in  the 
pitch  case,  the  resulting  equations  have  the  form 


+  4R  9  +  (R  -  1) 
x  x 


(81) 
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de2 
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(82) 


It  is  seen  that  no  first-order  terms  involving  the  eccentricity 
are.  present  and  thus,  no  yaw-roll  perturbation  is  caused  by  eccentricity. 
Rotating  Equipment 

If  rotating  equipment  is  included  in  the  payload  of  a  satellite, 
the  vehicle  can  no  longer  be  treated  as  a  single  rigid  body,  but  must 
be  considered  as  a  system  of  coupled  bodies,  consisting  of  the  main 
body  of  the  vehicle  and  the  various  internal  moving  parts.  Since 
angular  momentum  must  be  conserved,  it  is  obvious  that  the  motion  of 
these  internal  parts  can  react  on  the  main  satellite  body  to  produce 
disturbances  in  its  attitude. 

This  analysis  is  restricted  to  the  case  of  a  single  rotating  ele¬ 
ment  on-board  the  satellite.  It  is  assumed  that  the  element  has  ro¬ 
tational  symmetry  about  its  spin  axis  and  that  this  axis  lies  along 
one  of  the  principal  axes  of  the  main  satellite  body.  Under  these 
conditions,  the  motion  of  the  main  body  is  specified  by  Eqs.  (17) -(19) 

where  M  ,  M  and  M  represent  the  moments  applied  to  the  main  body  by 
x  y  z 

the  rotating  element  either  through  its  bearings  or  through  its  motor 


drive . 
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Pi.twh  Excitation.  If  the  rotating  element  has  i 
the  pitch  axis  of  the  main  vehicle,  its  equations  of 

by  d-,, 

h  dT  +  h  -  tt  “A  =  3“o<  V  V  ° 


ts  spin  axis  along 
motion  are  given 


(83) 
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s  z2x  o  T  s  vv 
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-  M 


z 


(84) 

(85) 


where  I  is  the  moment  of  inertia  of  the  rotating  element  about  its 
s 

spin  axis  and  1^  is  the  value  for  the  two  transverse  axes.  The 
quantity  cd^  is  the  total  angular  velocity  of  the  rotor  about  its  spin 
axis  and  is  given  by 


U) 


z2 


w  +  a. 


(86) 


where  is  the  spin  rate  of  the  rotor  relative  to  the  main  body. 

Substitution  of  Eqs „  (10)-(15)  and  Eq.  (86)  into  Eqs,  (83)-(85)  yields 
the  following  small  angle  equations  for  the  motion  of  the  rotor 

ITe  +  4w2(Ig-  IT)  9  +  <uo*(Is-  2It)  +  Isd2(#  +  aoB)  =  -  Mx  (87) 


IT^'  +  “o<V  *T^  ‘  fUo®(V  2IT^  ‘  1s*2(‘?  "  “  "  My 


(88) 


V®  +  »2)  =  "  Mz  (89) 

Since  there  is  no  differential,  rotation  of  the  two  bodies  about  the 

x  or  y  axes,  the  quantities  M  and  M  represent  constraint  moments 

x  y 

applied  through  the  bearings  normal  to  the  axis  of  rotation.  Thus, 
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M  and  M  can  be  eliminated  bv  adding  Eqs.  (17)  and  (87)  as  well  as 
x  y 

Eqs  ,  (18)  and  (88)  to  give 
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(91) 


(92) 


An  examination  of  Eqs.  (90)  and  (91)  shows  that  the  effect  of  is 
to  increase  the  yaw-roll  coupling  as  well  as  the  yaw  and  roll  stiffness. 
However,  if  the  initial  conditions  for  the  yaw  and  roll  angles  and  8 
are  zero,  the  angular  rate  produces  no  disturbance  in  either  yaw 
or  roll. 

Thus,  the  effect  of  ty^  Is  to  cause  variations  in  the  pitch  attitude 

of  the  main  body  as  described  by  Eqs.  (19)  and  (89).  In  this  case  the 

moment  M  is  exerted  through  the  motor  drive  on  the  two  bodies.  In 
z 

order  to  solve  Eqs.  (19)  and  (89),  it  is  necessary  to  obtain  an  ex¬ 
pression  for  as  a  function  of  at 2'  This  is  done  with  a  simple  motor 
equation  of  the  form 

-  M  =  k  i  -  for 
z  c 


(93) 
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where  -  M  ,  the  moment  applied  to  the  rotor,  is  made  up  of  a  terra  pro- 

Z  r 

portional  to  the  current  i,  and  a  viscous  drag  term  proportional  to  o,,. 
I^ie  motor-  current  can  be  expressed  by  the  relation 


E  =  R:  + 


(94) 


where  the  applied  voltage,  E,  is  equal  to  the  drop  across  the  motor 
resistance  plus  a  back  emf  term  proportional  to 

Substitution  of  Eqs.  (93)  and  (94)  into  Eqs .  (19)  and  (89)  gives 
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where 
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(97) 


It  can  be  shown  that  if  I  »  I  ,  then  the  solution  for  Eqs.  (95) 
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and  (96)  for  a  seep  change  in  motor  voltage  is  given  b> 
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where 
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(100) 
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and  is  the  time  required  for  the  motor  to  come  up  to  its  steady  speed. 
The  value  of  this  steady  speed,  ,  is  determined  from  Eq .  (99)  as 
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s 


and  Eqs,  (98)  and  (99)  can  be  put  in  the  form 
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If  the  time  constant  is  small  compared  to  the  orbital  period,  then 


u)  T  «■  1 
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and  Eq.  (102)  can  be  written  as 
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Thus,  if  the  voltage  applied  to  the  motor  is  increased  by  a  step, 
an  oscillation  in  the  pitch  attitude  of  the  main  body  will  result,  as 
indicated  in  Eq .  (104). 

Yaw  Excitation.  If  the  rotating  element  has  its  spin  axis  along 
the  yaw  axis  of  the  main  vehicle,  its  equations  of  motion  are  given 
by 
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ty 

IT  "rr — •  +  I  U)  ~u)  -  I_  tii  a)  =  3 <  (I  -  I_,)  a  b  -  M  (107) 

I  at  sy2x  Txy  os  T  yy  z  v  1 


As  before,  I  and  I  are  the  spin  and  transverse  moments  of  inertia 
s  T 

respectively,  and  tu  _  is  the  total  angular  velocity  of  the  rotor  about 

y  *- 

its  spin  axis  given  by 


(108) 


Substitution  of  Eqs .  (10) -(15)  and  Eq.  (108)  in  Eqs.  (1Q5)-(107)  gives 
the  small  angle  form  as 
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where  it  is  assumed  that  a  «  a>  ,  In  this  case,  M  and  M  are  constraint 

o  X  z 

moments  and  can  be  eliminated  between  Eqs.  (17)  and  (109)  and  Eqs.  (19) 
and  (111),  respectively,  to  give 
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An  examination  of  Eq.  (113)  shows  that  the  terra  involving  is 
negligible.  Thus,  the  pitch  behavior  of  the  vehicle  is  essentially 
independent  of  internal  rotations  about  the  yaw  axis, 

It  is  now  necessary  to  obtain  an  expression  for  the  interaction 
moment  M  .  This  is  done  in  the  same  manner  as  in  the  pitch  case  with 

y 

the  result  that 


M  = 

y 


(115) 


The  yaw-roll  behavior  of  the  vehicle  is  now  specified  by  Eqs.  (18), 
(110),  (111)  and  (115),  which  can  be  combined  to  give  the  following 
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(118) 


where  the  definitions  of  tu„  ,  u),  ,  and  p  are  of  the  same  form  as 

6  i!r  P  V 

Eqs.  (35)-(38),  using  the  primed  values  of  the  moments  of  inertia  to 

determine  R  and  R  . 

x  y 

It  can  be  shown  that  if  I  is  small  compared  to  the  orbital  period, 

c 

then  ti), T  and  tu_T  are  small  compared  to  unity  since  from  the  Appendix 
1  o  2  o 

it  is  seen  that  and  are  of  the  same  order  as  o>o.  If.,  in  addition, 
I  is  small  compared  to  both  i'  and  I7,  then  a  good  approximation  to 
the  solution  of  Eqs.  (116) -(118)  for  a  step  change  in  E  is  given  by 
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where  the  steady-state  angular  rate  of  the  motor,  |0  ,  is  given  >y 
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and 
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(123) 


The  values  of  and  ^  are  determined  as  described  in  the  Appendix, 
using  the  primed  values  cf  the  moments  of  inertia  to  determine  and 

R  . 

y 

It  is  seen  that  the  exponential  terms  in  Eqs .  (119) -(121)  become 
negligible  after  a  relatively  short  time  (of  the  order  of  a  few  time 
constants).  Thus,  the  roll  response,  Eq ,  (120),  is  made  up  of  a 
steady-state  term  and  two  oscillatory  terms  of  frequencies  ^  and  ^ 
The  amplitude  of  the  steady-state  terr  is  given  by 
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The  oscillatory  terms  in  the  rv  -  response  have  a  maximum  possible 
amplitude  equal  to  the  sum  of  their  individual  amplitudes  given  by 
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This  amplitude  is  equal  to  the  steady-state  value  determined  in 
Eq .  (124).  Thus,  the  maximum  roll  angle  would  be  twice  that  indicated 
in  Eqs .  (124)  c-r  (125).  However,  in  the  presence  of  damping  the  os¬ 
cillatory  amplitude  would  reduce  to  zero  while  the  steady-state  value 
would  remain  unaffected. 

The  yaw  response  consists  of  two  oscillatory  terms  with  a  maximum 
possible  amplitude  given  by 
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(126) 


Again,  it  should  be  emphasized  that  a  damping  system  would  reduce  this 
yaw  oscillation  to  zero. 
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Roll  Excitation.  When  the  rotating  element  has  its  spin  axis  along 
the  roll  axis  of  the  main  vehicle,  the  equations  of  motion  are  developed 
in  a  completely  analogous  manner  to  that  used  in  the  case  of  yaw  exci¬ 
tation.  The  resulting  pitch  equation  has  the  form 
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As  is  the  yaw  case  above,  the  term  in  Eq.  (127)  involving  the  angular 
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rate  of  the  motor,  Is  negligible,  and  the  pitch  equation  is  essen¬ 
tially  independent  of  rotations  about  the  roll  axis. 

The  equations  describing  the  motion  of  the  vehicle,  in  yaw  and  roll 
and  the  coupled  motion  of  the  motor  have  the  form 
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under  the  assume tion  that  y  «  w  . 

“  o 

If,  as  before,  (1,1  and  u.„T  are  much  less  than  unity  and  I  is 
>  7  o  L  o  s 

much  less  than  l'  and  l'  ,  then  the  solution  to  Eqs.  (129) -(131)  for 
x  y 


a  step  in  E  is  given  by 
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(1-34) 
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As  before,  the  exponential  terms  in  Eqs.  (133)  and  (134)  can  be  neg¬ 
lected  after  a  short  time,  and  it  is  seen  that  the  steady-state  yaw 
response  is  given  by 
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The  maximum  possible  amplitude  of  the  oscillatory  part  of  the  yaw  re¬ 
sponse  is  the  sum  of  the  absolute  values  of  the  amplitudes  of  the  last 
two  terms  of  Eq .  (134)  expressed  in  the  form 
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In  this  case  it  is  seen  that  is  greater  than  s s  - 

The  maximum  possible  amplitude  of  the  oscillatory  roll  response 
is  given  by 
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III.  RESULTS 

In  this  section,  some  typical  values  of  the  attitude  perturbations 
of  a  satellite  vehicle  are  determined  by  means  of  the  formulas  developed 
in  the  p rev i  .s  section. 

ASSUMED  VEHICLE  CONFIGURATION 

While,  it  is  not  the  purpose  of  this  Memorandum  to  design  a  satellite 
vehicle,  it  is  necessary  to  specify  the  three  principal  moments  of 
inertia  of  the  vehicle  before  the  formulas  of  the  previous  section 
can  be  evaluated.  Thus ,  for  a  nominal  design  the  moments  of  inertia 
are  specified  as 

I  -  9500  slug  ft2 
x 

I  =  1000  slug  ft2  (140) 

I  =  10000  slug  ft2 

These  values  can  be  achieved  by  a  suitable  distribution  of  a  total 
vehicle  mass  of  750  lbs.  This  choice  of  moments  of  inertia  satis¬ 
fies  the  physical  restrictions  as  given  by  Eqs.  (186)-(188)  in  the 
Appendix. 

From  Eqs.  (20),  (21)  and  (22),  the  values  of  R^,  and  R^  are 
obtained  as 

R  =  .9474 
x 

R  -  .5 

y 

R  =  .85 
z 

It  is  seen  that  these  values  also  satisfy  the  conditions  for  pitch, 
yaw  and  roll  stability  as  specified  by  Eqs.  (192),  (195)  and  (196). 


30 


By  means  of  Eqr. .  (191),  the  characteristic  frequency  in  pitch  is 
found  to  be 


Oj  = 
O' 

1.5969  ui 

o 

(142) 

while  Eqs.  (193)  and  (194)  give 

the  characteristic  yaw-roll 

frequencies 

as 

\  = 

1.9531 

0 

(143) 

a) 2  = 

.7043  o) 

o 

(144) 

EFFECT  OF  AN  IMPULSE  MOMENT 

The  expressions  for  the  attitude  perturbations  due  to  an  impulse 
moment  were  determined  in  Section  II  and  will  be  applied  here  to  the 
case  of  the  impulse  resulting  from  meteoroid  impact. 

Pitch  Impulse 

The  maximum  pitch  attitude  perturbation  due  to  an  impulse  moment 
is  given  by  Eq .  (30),  which  when  combined  with  Eq .  (16)  gives 


a 

m 


.13225 


R1/2  X 


where  the  units  are  as  follows 


r  -  n  mi 
c 


I  -  ib-ft-sec 
P 

Iz  -  slug  ft^ 


cy  -  deg 
m 


(145) 


If  this  result  is  applied  to  the  case  of  meteoroid  impact,  the 


impulse  can  be  expressed  in  the  form 
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I  =  r  A  (mV)  (146) 

P  P 

where  A(mV)  the.  change  of  momentum  of  the  meteoroid,  which  is  equal 

to  the  force  impulse  in  lb-sec,  while  r  is  the  pitch  moment  arm  of 

P 

this  force  in  ft. 

The  determination  of  the  momentum  change  involved  in  meteoroid 
impact  is  extremely  complicated  and  involves  not  only  the  mass  and 
velocity  of  the  incident  particle  but  also  the  mass  and  velocity  of 
the  vaporized  target  material.  For  the  purposes  of  this  Memorandum, 
the  momentum  change  is  evaluated  as  the  product  of  the  mass  and  velocity 
of  the  incident  particle.  This  assumes  an  inelastic  impact  and  gives 
a  lower  limit  for  the  value  of  A(mV) . 

As  an  example,  suppose  the  satellite  vehicle  is  at  synchronous 
altitude  (r  =  2275Q  n  mi)  and  is  subjected  to  a  meteoroid  impact  for 
which 

'  ;  .,V)  =  10 "3  lb -sec 

r  =  10  ft 
P 

The  resulting  pitch  attitude  perturbation  has  a  maximum  value  given 
by  Eqs.  (145)  and  (146)  as 

ex  =  ,492  deg 
m 

It  should  be  noted  that  a  meteoroid  impact  corresponding  to  as 

_3 

much  as  10  lb-sec  would  occur  relatively  infrequently. 

Roll  Impulse 

The  maximum  roll  attitude  perturbation  due  to  an  impulse  moment 
in  roll  is  given  by  Eq.  (41)  which  can  be  expressed  as 
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and  the  units  are  the  same  as  those  used  in  the  pitch  case  above. 

As  indicated  in  Section  III,  the  coupling  between  yaw  and  roll 
results  in  a  yaw  response  even  though  I  is  applied  ai>out  the  roll 
axis.  The  amplitude  of  this  yaw  response  is  given  by  £q.  (42),  which 
can  be  written  in  the  form 
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.4580 


f2(V 


where 


F2<V  ®2>  ‘  ' 


(R  -  l) 

-2 . . 


049) 


(150) 


and 

IR  =  rRA^nV) 


(151) 


As  an  example,  suppose  that  a  meteoroid  impact  occurs  on  the  same 
synchronous  satellite  and  produces  an  impulse  moment  in  roll.  If  I 

R 

is  specified  by 

A (mV)  =  10'3  lb-sec 


rR  -  10 
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where  r  is  the  roil  moment  arm,  then  the  resulting  amplitudes  of  the 
R 

roll  and  yaw  perturbations  as  determined  by  Eqs.  (147)  through  (130) 
are  given  by 

j5  =  .425  deg 

m 

A  =  .249  deg 

•  m 


Yaw  Impulse 

The  maximum  yaw  attitude  perturbation  due  to  an  Impulse  moment  in 
yaw  is  given  by  Eq ,  (48),  which  can  be  expressed  as 
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The  amplitude  of  the  associated  roll  response  is  given  by  Eq.  (49), 
which  can  be  expressed  as 
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and 

ly  =  r^CmV)  (156) 

As  an  example,  if  an  Impulse  moment  due  to  a  meteoroid  impact  is 
characterized  by 

A  (mV)  =  10  ^  lb -sec 
ry  =  10  ft 

the  resulting  yaw  and  roll  amplitudes  as  determined  by  Eqs.  (152) 
through  (156)  are 

i t  =  11.11  deg 
m 

p  =  .249  deg 

m 

A  comparison  of  die  results  for  roll  and  yaw  impulses  shows  that  the 

largest  response  is  that  in  yaw  due  to  a  yaw  impulse.  It  is  also  seen 

that  g  due  to  a  yaw  impulse  is  equal  to  4  due  to  a  roll  impulse, 
m  Tm 

This  can  be  demonstrated  by  a  comparison  of  Eqs.  (42)  and  (49). 

EFFECT  OF  A  STEP  MOMENT 

The  expressions  for  the  attitude  perturbations  due  to  a  step  in¬ 
crease  in  the  applied  moment  were  determined  in  Section  II,  and  will 
be  applied  here  to  the  case  of  the  step  change  in  the  radiation  pres¬ 
sure  moment  which  occurs  as  the  vehicle  emerges  from  the  earth's  shadow. 
Step  Moment  in  Pitch 

The  steady-state  deflection  in  pitch  due  to  a  step  moment  in  pitch 
is  given  by  Eq .  (54),  which  can  be  rewritten  as 

-4  r3  Mp 

'  3'051  *  10 

Z  Z 


(157) 
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where  the  units  ai*c  as  follows 

r 

c 

I 

z 

«p 
n 

ss 

To  determine  the  magnitude  of  Mp,  assume  the  radiation  is  incident 

2 

along  the  x  axis,  and  that  the  satellite  presents  an  area  of  A  ft  in 

the  x  direction.  If  the  surface  is  totally  reflecting,  the  total  force 

-7  2 

due  to  radiation  pressure  (p  =  1.865  x  10  lbs/ft  )  is 

FR  ’  ”0*  <158> 

acting  along  the  x  axis.  If  the  center  of  pressure  is  displaced  a 
distance  rp  up  the  y  axis,  the  resulting  pitch  moment  is  expressed  as 

Mp  -  PQ  Arp  (159) 

7 

As  an  example,  if  the  subtended  area  is  100  ft”  and  the  center  of 
pressure-center  of  mass  offset,  Cp,  is  .1  ft,  than  the  pitch  moment 
is  given  by 

Mp  =  1  .865  x  10"b  lb-ft 

If  a  step  change  in  moment  of  this  magnitude  were  applied  to  the  assumed 
satellite  configuration  in  a  synchronous  orbit,  the  resulting  steady- 
state  pitch  angle  based  on  Eq .  (157)  would  be 

a 

ss 


n  mi 

2 

slug  ft“ 

lb-ft 

deg 


=  .788  deg 


36 


Step  Moment  in  Roll 

The  steady -state  deflection  in  roll  due  to  a  step  moment  in  roll 
is  g’ven  by  Eq .  (61),  which  can  be  written  in  the  form 


ss 


=  2.288  x  10 
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In  addition,  the  amplitude  of  the  a-  eclated  yaw  response  is  given 
by  Eq .  (60),  which  can  be  expressed  as 


where 
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r3  M 

=  9.153  x  10  G  (t»  ,  -x2) 
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(161) 


(162) 


If  the  step  moment  Mp  has  a  magnitude  equal  to  that  used  for  Mp, 
then  the  resulting  response  of  the  satellite  as  described  by  Eqs.  ( 160) — 
(162)  is 

8  -  .558  dee 

^  =  .616  deg 

m 


Step  Moment  in  Yaw 

The  steady-state  deflection  in  yaw  due  to  a  step  moment  in  yaw 
is  given  by  Eq .  (68),  which  can  be  rewritten  in  the  form 


$  =  9.153  x  10"4 

ss 
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The  amplitude  of  the  associated  roll  oscillation  is  given  by 
Eq .  (67),  which  can  be  expressed  as 
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If  the  amplitude  of  the  step  moment  in  yaw  is  the  same  as  that 
assumed  ^or  the  other  two  axes,  Eqs.  (163)  and  (164)  give  the  follow¬ 
ing  results 

iss  =  40,2  deg 

8  =  .616  deg 

m 

Again,  it  is  seen  that  the  system  is  extremely  sensitive  in  its  yaw 
response  and  exceeds  the  small  angle  assumption  used  in  the  derivation 
of  the  equations  of  motion.  In  addition,  the  yaw  amplitude  resulting 
from  a  step  moment  in  roll  is  equal  to  the  roll  amplitude  resu1 ting 
from  a  step  moment  in  yaw.  This  can  be  shown  by  a  comparison  of 
Eqs.  (cO)  and  (6>). 

In  general,  these  results  emphasize  the  importance  of  maintaining 
the  certer  ol  radiation  pressure  of  the  vehicle  close  to  the  center  of 
mass . 

EFFECT  OF  STATION-KEEPING  PROPULSION 


In  many  satellite  applications  it  may  be  necessary  to  make  small 
adiustments  in  velocity  in  order  to  maintain  station  relative  to  other 
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satellites  or,  as  in  the  case  of  a  synchronous  satellite,  to  keep  the 
same  position  relative  to  the  earth.  If  the  thrust  used  to  produce 
these  velocity  changes  has  a  line  of  action  which  passes  through  the 
vehicle  center  mass,  r>o  attitude  disturbance  will  result.  However, 
determination  of  the  permissible  misalignment,  is  an  important  design 
problem. 

Impulsive  Velocity  Change 

It  is  assumed  that  station-keeping  adjustments  are  restricted  to 
velocity  changes  along  the  orbit.  Thus,  the  resulting  moments  are 
about  the  pitch  or  yaw  axes.  If  the  velocity  adjustment  is  made  im¬ 
pulsively,  then  the  amplitude  of  the  resulting  pitch  oscillation  can 
be  determined  by  Eqs.  (145)  and  (146).  However,  in  this  case, 
the  quantity  A (mV)  represents  the  desired  momentum  change  of  the  ve¬ 
hicle  and  Tp  is  a  measure  of  the  thrust  misalignment.  Similarly,  the 
amplitude  of  the  yaw  response  and  ics  associated  roll  response  are 
given  by  tqs.  (152)-(156)  with  the  above  definition  of  £(mV). 

As  an  example,  suppose  th*.t  the  assumed  satellite  vehicle,  weighing 
750  lbs,  is  in  a  synchronous  orbit  with  a  radius  of  22,750  n  mi.  If  a 

velocity  change  of  I  ft/sec  is  required  and  the  moment  arm,  r^,  of  the 

-3 

thrust  about  the  pitch  axis  is  equal  to  10  ft,  then  the  amplitude  of 
the  resulting  pitch  oscillation  as  given  by  Eq .  (145)  is 

a  -  1 . 146  deg 
tn 

Similarly,  if  the  thrust  has  a  moment  arm,  r^,  about  the  yaw  axis  equal 
Lo  10  ft,  then  the  amplitudes  of  the  resulting  yaw  and  roll  oscillations 
as  given  by  Eqs.  (152)  and  (134)  are 
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OS  a-r 

W  t 


deg 


y  =  .58  deg 

m 


The  assumed  moment  arm  of  10  ft  represents  a  precision  of  thrust 
alignment  which  is  not  achievable  in  practice. 

Constant  Thrust  Application 

In  Ref.  4  the  station-keeping  propulsion  requirements  resulting 
from  the  ellipticity  of  the  earth's  equator  was  determined  as  17  ft/ 
sec/yr.  However,  based  on  observations  of  Sy.com  2,  a  more  realistic 

value  appears  to  be  6  ft/sec/yr.  For  the  assumed  750-lb  satellite 

-6 

this  requirement  is  equivalent  to  a  steady  thrust  of  5.17  x  10  lb. 

If  the  moment  arm,  rp,  of  this  force  about  the  pitch  axis  is  .1  ft, 
then  the  resulting  steady-state  pitch  angle  can  be  determined  by  means 
of  Eq.  (157)  as 


a  =  .218  deg 

ss 

Similarly,  if  the  moment  arm,  r^,  of  this  force  about  the  yaw  axis  is 
.1  ft,  then  the  resulting  steady-state  yaw  angle  as  given  by  Eq .  (163) 
is 

if  =11.14  deg 

'  SS 

Thus,  it  is  seen  that  even  at  extremely  low  thrust  levels,  appreciable 
attitude  disturbances  can  occur. 

EFFECT  OF  ORBITAL  ECCENTRICITY 

The  amplitude  of  the  forced  oscillation  in  pitch  due  to  orbital 
eccentricity  is  given  by  Eq .  (78),  which  can  be  expressed  as 
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where  ry  is  in  degrees, 
m 

As  an  example,  if  the  assumed  satellite  configuration  is  or  an  or¬ 
bit  with  an  eccentricity  of  .005,  Eq .  (166)  gives  a  pitch  amplitude  of 


0!  =  .370  deg 

m 

An  examination  of  Eq .  (78)  shows  that  the  value  of  jy  depends  only  on  the 
satellite  configuration  and  the  orbital  eccentric it,  and  is  independent 
of  the  orbital  period.  It  is  seen  that  <y  in  radians  is  approximately 
equal  to  the  eccentricity;  thus  a  low  eccentricity  orbit  is  essential 
in  order  to  avoid  large  errors  in  pitch  attitude. 

EFFECT  OF  ROTATING  EQUIPMENT 

In  Section  II,  the  expressions  were  developed  for  the  attitude  per¬ 
turbations  which  would  result  due  to  the  presence  of  rotating  equipment 
on-board  a  satellite.  These  relations  are  used  here  to  determine  repre¬ 
sentative  values  of  these  perturbations  for  the  assumed  satellite  con¬ 
figuration. 

Rotation  about  Pitch  Axis 

From  Eq .  (104)  it  is  seen  that  if  a  rotor  of  moment  of  inertia  I 

s 

is  accelerated  to  a  steady  rotation  rate  of  at-  relative  to  the  main 

^  s  s 

vehicle,  an  oscillation  in  pitch  is  initiated  which  has  an  amplitude 
given  by 


(i67) 
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where  Gp  is  the  angular  momentum  ot  the  rotor. 


Equation  (167)  can  be 


rewritten  as 


3/2 
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i. 

a  =  .13221  -|t t 

m  AH 
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(168) 


where  the  units  are 


n  mi 


I  -  slug  ft' 


a 


m 


slug  ft  /sec 
deg 


The  magnitude  of  Gp  is  given  by  the  relation 


G  =  .05236  I  d„ 

P  ?  2s  s 


(169) 


where  the  units  are 


slug  ft 


ot. 


rpm 


slug  ft  /sec 


As  an  example,  suppose  the  assumed  satellite  configuration  is  on  a 

2 

synchronous  orbit  and  has  a  rotor  with  1^  =  .001  slug  ft  and  a  rotation 
rate  ~  100  rPm 

By  means  of  Eqs.  ( 16 S)  and  (169)  the  amplitude  of  the  resulting  pitch 
oscillation  is  found  to  be 


oq  =  .258  deg 
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It  should  be  noted  that  if  an  adequate  damping  system  is  present, 
this  pitch  oscillation  will  be  reduced  to  zero.  However  each  time 
the  rotor  starts  or  stops  a  similar  transient  oscillation  will  be 
generated . 

Rotation  about  Yaw  Axis 

If  the  internal  rotation  is  tbout  the  yaw  axis,  the  resulting  yaw 
response  is  oscillatory  with  an  amplitude  determined  by  Eq .  (126), 
which  can  be  written  in  the  form 


in 


.22904 
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where 
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In  addition,  the  associated  roll  response  consists  of  a  steady- 
state  deflection  plus  an  oscillatory  term  with  an  amplitude  equal  to 
the  steady-state  term.  The  magnitude  of  the  steady-state  term  is  given 
by  Eq .  (124),  which  can  be  written  as 


a 


ss 


.05726 


(173) 
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If  the  value  ot  G^,  is  identical  to  that  chosen  for  ,  the  re¬ 
sulting  yaw  and  roll  responses  of  the  assumed  synchronous  satellite 
are.  given  by 

I?  =  5.876  deg 
fm 

6  =  .114  deg 

m 

g  =  .114  deg 

ss 

As  in  the  pitch  case,  suitable  damping  will  reduce  the  oscillatory 
amplitudes  to  zero.  However,  the  steady-state  roll  angle  is  unaffected 
as  long  as  the  internal  rotation  continues. 

Rotation  about  Roll  Axis 

If  the  internal  rotation  is  about  the  roll  axis,  the  resulting  roll 
response  is  oscillatory  with  an  anplitude  given  by  Eq.  (139),  which 
can  be  written-  in  the  form 


3 


m 


,22904 


where 


1^2  C  a)  i  >  u^) 


tiJl 

i  +  -M- 

t. 


and 


GR  “  *°5236  hss 


(174) 


(175) 


(176) 


In  audition,  the  associated  yaw  response  is  made  up  of  a  steady- 
state  term  given  by  Eq.  (137),  which  can  be  expressed  as 
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*ss 


.22904 


(1/7) 


and  also  an  oscillatory  term  with  an  amplitude  given  by  Eq  .  (138),  which 
can  be  put  in  the  form 


3/2 


ijr  = 

m 


.22904 


i'- 


1' 

x 


tii/- 


(178) 


where 


H3G01,  ii  2 )  " 


(179) 


If  the  value  cf  G  is  the  same  as  that  assumed  for  G  and  G.r,  the  yaw- 

K.  r  i 

roll  response  of  the  satellite  as  given  by  Eqs.  (174),  (177)  and  (178) 
is  as  follows 

B  =  .282  deg 

m 

$  -  8 . 34  deg 

«ss  =8.23  cleg 

As  in  the  previous  cases,  the  presence  of  damping  will  reduce  the 
oscillatory  amplitudes  to  zero  after  a  few  orbits,  but  the  steady-state 
yaw  angle  will  remain  as  long  as  the  rotor  maintains  its  steady-state 
rate  e2ss. 
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IV.  DISCUSSION  AND  CONCLUSIONS 

An  examination  of  the  results  presented  In  Section  III  indicates 
that  a  satellite  which  depends  on  gravity  gradient  to  maintain  attitude 
stability  may  be  subject  to  appreciable  perturbations.  This  results 
from  the  fact  that  although  the  principal  moments  of  inertia  can  be 
chosen  so  that  the  vehicle  is  inherently  stable,  the  restoring  torques 
are  extremely  small.  Thus ,  it  is  necessary  to  consider  the  effects  of 
very  small  disturbances  such  as  meteoroid  impact,  radiation  pressure, 
station-keeping  propulsion,  orbital  eccentricity  and  internal  rotating 
equipment  which  would  ordinarily  be  ignored  ir  the  presence  of  an 
active  attitude  control  system. 

The  results  indicate  that  a  meteoroid  impact  capable  of  producing 
an  impulsive  moment  equal  to  .01  lb-ft-sec  can  produce  deflections  of 
a  synchronous  satellite  of  the  order  of  a  fraction  of  a  degree  in  either 
pitch  or  roll  and  of  several  degrees  in  yaw.  Hie  frequency  with  which 
such  an  impact  might  occur  depends  on  the  mechanism  assumed  for  hyper¬ 
velocity  impact  and  a  more  precise  description  of  the  vehicle  con¬ 
figuration.  Thus,  the  results  are  given  in  terms  of  the  momentum 
change  occurring  at  impact. 

In  the  case  of  radiation  pressure  it  is  seen  that  for  the  assumed 
satellite  configuration  at  synchronous  altitude,  a  moment  of  1.865  x  10 
Ib-ft  due  to  radiation  pressure  could  cause  deflection  of  the  order  of 

a  fraction  of  a  degree  in  pitch  and  roll  and  tens  of  degrees  in  yaw. 

2 

The  assumed  moment  corresponds  to  the  force  on  100  rt  of  reflecting  sur¬ 
face  with  a  moment  arm  of  .1  ft.  This  emphasizes  the  importance  of 
maintaining  symmetry  of  the  reflecting  properties  of  the  vehicle. 
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Impulsive  station-keeping  propulsion  requires  extreme  precision  in 
the  alignment  of  the  thrust  axis  through  the  center  of  mass.  Even  for 
extremely  low  thrust  levels  and  long  burning  times  the  precision  is 
still  excessive.  Thus,  it  appears  inadvisable  to  use  a  gravity  gradient 
system  if  s tat ion -keeping  is  required. 

Orbital  eccentricity  causes  an  oscillation  in  pitch  of  .37  deg  for 
an  eccen  tricity  of  .005.  This  eccentricity  corresponds  to  a  difference 
of  about  225  n  mi  in  apogee  and  perigee  altitudes  for  a  synchronous  or¬ 
bit.  Actually,  Syncom  and  Early  Bird  have  achieved  appreciably  lower 
eccentricities . 

In  the  case  of  rotating  equipment,  it  is  seen  that  an  internal 

piece  of  rotating  equipment  with  an  angular  momentum  of  .005236  slug- 
2 

ft  /sec  aligned  along  the  pitch  axis  produces  an  oscillatory  deflection 
in  pitch  of  about  .25  deg.  The  same  equipment  aligned  along  the  yaw 
axis  produces  a  steady-state  roll  angle  of  about  .11  deg  and  a  yaw 
oscillation  of  5.88  deg.  Finally,  if  the  axis  of  rotation  is  along  the 
roll  axis  a  steady-state  yaw  angle  of  8.23  deg  is  produced  as  well  as 
a  roll  oscillation  of  about  .28  deg. 

To  summarize  the  results  presented  above,  the  various  sources  of 
perturbation  are  listed  below  in  the  order  of  their  importance  and  with 
indications  of  how  their  effects  can  be  reduced. 

1 .  Station -Keeping  Propulsion 

It  is  extremely  doubtful  that  the  required  precision  of  thrust 
alignment  can  be  achieved  in  practice.  Thus,  Station-Keeping  is  in¬ 
compatible  with  gravity  gradient  stabilization,  unlesj  auxiliary  sen¬ 


sors  are  used  to  determine  the  attitude. 


47 


2.  Rotating  Equipment 

If  such  equipment  is  essential  to  the  success  of  the  mission,  its 
perturbing  effect  can  be  reduced  by  utilizing  similar  counter-rotating 
elements.  Thus ,  the  quantities  G^ ,  G^,  and  are  the  residual  un¬ 
compensated  angular  momenta. 

3 .  Radiation  Pressure 

In  order  to  minimize  the  effects  of  radiation  pressure,  it  is 
necessary  to  maintain  geometrical  symmetry  as  well  as  symmetry  of  the 
surface  reflectivity  so  that  the  force  due  to  radiation  acts  through 
the  center  of  mass  of  the  vehicle.  Hiis  may  be  difficult  to  achie\e 
if  the  vehicLe  is  subject  to  thermal  bending  and  resultant  variations 
in  the  position  of  the  center  of  mass. 

4 .  Meteoroid  Impact 

The  relative  importance  of  the  effect  depends  upon  the  statistical 
distribution  of  meteoroid  impacts  as  a  function  of  impulse  as  well  as  a 
detailed  description  of  the  vehicle  configuration,  both  of  which  are 
beyond  the  scope  of  this  Memorandum. 

5 .  Orbital  Eccentricity 

Th  .s  effect  requires  that  the  vehicle  be  on  a  nearly  circular  orbit 
to  avoid  a  steady  oscillation  in  pitch.  However,  eccentricities  cur¬ 
rently  being  achieved  are  lov;  enough  that  this  effect  should  not 
be  serious. 

Finally,  it  should  be  pointed  out  that  for  many  applications  attitude 
variations  of  the  vehicle  of  the  order  of  a  few  degrees  may  be  acceptable 
if  these  variations  can  be  measured  precisely  and  read  out  to  the  user. 
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Append  ix 

DETERMINATION  OF  VEHICLE  CONFIGURATION  AND  STABILITY 

In  Section  II  of  this  Memorandum,  the  linearized  equations  of  motion 

were  developed  in  terms  of  the  parameters  R  and  R  (see  Eqs.  (23)-(26)), 

x  y 

and  it  was  pointed  out  that  there  are  certain  restrictions  on  the  per¬ 
missible  values  of  R  and  R  . 

x  y 

PHYSICAL  RESTRICTION 

The  fact  that  the  three  principal  moments  of  inertia  cannot  be 
selected  arbitrarily  puts  a  limitation  on  the  values  of  R  and  R 

x  y 

through  Eqs.  (20)  and  (21)  as  follows. 

The  three  moments  of  inertia  can  be  expressed  in  the  form 
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f  y2dm  +  [  z 
n  JM 


z^dm 
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[  y2dm 
JM 


(182) 


If  these  three  equations  are  solved  simultaneously  for  the  three 
integrals,  their  values  are  given  by 


[  x2dm  = 
JM 


I  +  I  +  I 
x  y  z 
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2  r  Z2dm  =  I  +  I  -  I  (185) 

x  V  2 

Since  the  left  sides  of  these  equations  are  always  positive,  the 
following  conditions  must  be  satisfied  by  the  moments  of  inertia 


+  I  s  I 

(186) 

Z  X 

+  I  SI 

(187) 

x  y 

+  I  si 
y  z 

(188) 

Thus,  it  is  seen  that  the  sum  of  any  two  principal  mcments  of  inertia 
is  greater  than  or  ecual  to  the  third. 

By  means  of  the  above  inequalities  and  the  definitions  of  Rx  and 
Ry ,  Eqs.  (20)  and  (21),  it  can  be  shown  that 


-  1  <;  R  s  +  1 
x 


-  1  s  R  «;  +  1 

y 


(189) 

(190) 


Thus,  in  the  R  ,  R  plane,  all  physically  realizable  configurations 
x  y 

lie  inside  of  the  square  shown  in  Fig.  2. 


STABILITY  RESTRICTION 

In  addition  to  the  physical  restriction  specified  above,  it  is  also 

necessary  to  rule  out  combinations  of  R  and  R  for  which  the  vehicle 

x  y 

is  unstable  in  either  pitch,  yaw  or  roll. 


Pitch  Stability 

An  examination  of  Eq.  (25)  shows  that  the  pitch  response  is  com¬ 
pletely  uncoupled  from  the  yaw  and  roll  motion  and  has  a  characteristic 


frequency  given  by 
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R  R 
x  y 


(191) 


Since  the  absolute  values  of  R^  and  R^  are  less  than  unity,  Eq.  (191) 
represents  an  additional  restriction  such  that 


R  L*  R  (192) 

x  y  v  ' 

in  order  that  the  vehicle  be  stable  in  pitch.  This  condition  rules 

out  the  region  in  the  R  ,  R  plane  above  the  45°  line  (R  *  R  )  (see 

x  y  x  y 

Fig.  2). 

Yaw-Roll  Stability 

An  examination  of  Eqs.  (23)  and  (24)  shows  that  the  yaw  and  roll 
responses  are  counled  and  if  the  external  moments  and  are  set 
equal  to  zero,  the  two  natural  frequencies  for  the  yaw- roll  motion  can 
be  determined  from  the  characteristic  equation  as 


(193) 


“2 

T 

tU 


(3R  +  R  R  +  1) 
x  x  y 


(3R  +  R  R  +  1) 
x  x  y 


16  R  R 
_ XJL 


(194) 


These  two  equations  impose  additional  restrictions  on  R^  and  R^  since 

2  2 

it  is  necessary  that  and  6e  positive  and  real.  The  first  of  these 


conditions  is  satisfied  if  R  and  R  have  the  same  sign,  that  is 

x  y 
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R 
x  y 


;>  0 


(195) 


Ihis  condition  rules  out  the  second  and  fourth  quadrant  of  the  R  , 

X 

R  plane  although  the  second  quadrant  has  already  been  ruled  out  by 
Eq.  (192). 

2  2 

The  requirement  that  aj,  and  tu0  be  real  is  satisfied  if 

1  2 


(3R  +  R  R  +  l)2 
y  x  y 


^  16  R  R 
x  y 


(196) 


which  rules  out  the  region  in  the  third  quadrant  of  the  R^,  R^  plane 
below  the  curve 

(3R  +  R  R  +  l)2  =  16  R  R  (197) 

xxy  xv  x 

as  shown  in  Fig.  2. 

R  ,  R  PLANE 
x*  y 

The  physical  and  stability  restrictions  developed  above  have  de¬ 
fined  the  permissible  region  of  the  R^,  Ry  plane  as  shown  in  Fig.  2. 
This  is  similar  to  the  result  obtained  by  DeBra  and  Delp  in  Ref.  5. 

The  triangular  portion  of  the  region  in  the  first  quadrant  corre¬ 
sponds  to  those  configurations  ordinarily  considered  for  gravity  gradi¬ 
ent  stability  in  which 


I  >  I  >  I 
z  x  y 


(198) 


The  origin  (R^  =  Ry  =  0)  corresponds  to  a  spherical  mass  distri¬ 
bution  while  the  hypotenuse  (R  =  R  )  represents  configurations  which 

x  y 

have  rotational  symmetry  about  the  z  or  pitch  axis  (1^  *  1^).  The 
horizontal  edge  of  the  triangle  (Ry  *  0)  corresponds  to  configurations 
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with  rotational  symmetry  about  the  y  or  yaw  axis  (1=1).  Finally, 

2*  X 

the  vertical  edge  of  the  triangle  (R  =  1)  represents  any  linear  mass 

X 

distribution  along  the  y  or  yaw  axis  (I  =  0). 

It  should  be  noted  that  the  point  =  1  represents  a  planar 

mass  distribution  in  the  x,  y  plane  (I  +  I  =  I  ).  The  exact  nature 

x  y  z 

of  this  mass  distribution  may  be  anything  between  a  circular  disc  and 
a  linear  distribution  along  the  y  axis. 

The  other  portion  of  the  R^,  R^  plane  which  is  permissible  lies 
in  the  third  quadrant  as  shown  in  Fig.  2.  It  is  bounded  by  the  negative 
R  axis:  th'j  line,  R  =  R  :  the  line,  R  =  -  1 :  and  the  curve  specified 
by  Eq.  (197).  In  this  area  the  three  moments  of  inertia  have  the 
following  relation 


I  >  I  >  I 
x  y  z 


(199) 


with  I  =  I  along  the  45  line  (R  =  R  )  and  I  =  I  along  the  negative 
x  y  x  y  y  z 

Ry  axis  (Rx  =0).  It  is  seen  that  the  stability  in  this  case  is  that 
of  an  elongated  body  spinning  about  its  axis  of  minimum  moment  of  in¬ 
ertia  at  orbital  angular  rate. 

Pitch  Frequency  Contours 

It  is  of  interest  to  determine  the  contours  of  constant  pitch  fre¬ 
quency  on  the  R^,  plane.  This  can  be  done  by  means  of  Eq.  (191). 

The  resulting  contours  are  plotted  in  Fig.  3.  It  is  seen  that  in  both 
parts  of  the  st"  .. o  region  the  value  of  is  limited  to  the  range 


(JU 


0  *  JL  * 


U), 


(200) 


The  dotted  portions  of  contours  correspond  to  regions  of  yaw-roll  in¬ 
stability.  Thus,  although  the  pitch  motion  is  apparently  stable,  the 


Fig. 3— Contours  of  constant  pitch  frequency 
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linearize  1  equation  describing  the  pitch  motion,  Eq.  (25),  would  no 
longer  be  valid  in  the  presence  of  large  yaw  and  roll  angles. 
Yaw-Roll  Frequency  Contours 

The  characteristic  yaw-roll  frequencies  specified  by  Eqs.  (193) 
and  (194)  are  roots  of  the  equation 


%-  -  Jfijr  (3R  +  R  R  +  1)  +  4R  R  = 
4  2  x  xy  xy 


(201) 


This  function  can  be  used  to  determine  contours  of  constant  yaw-roll 

frequency  in  the  R^,  R  plane.  These  contours  are  plotted  in  Fig.  4 

x  y 

while  Fig.  5  is  an  enlargement  of  the  stable  region  in  the  third  quad¬ 
rant. 

An  examination  of  these  two  figures  shows  that  through  any  point 

of  the  stable  region  of  the  R  ,  R  plane  there  are  contours  corresponding 

x  y 

to  two  values  of  frequency.  These  two  values  are  the  and  tu  2  de‘ 
fined  by  Eqs.  (193)  and  (194).  It  is  seen  that  in  the  first  quadrant 


1  <  —  <  2 
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(202) 


0  <  —  <  1 


(203) 


while  in  the  third  quadrant  both  frequencies  are  less  than  u>  so  that 
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'°2  ^1  , 

—  ■'  —  <1 
0!  (D_ 
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(204) 


It  can  also  be  shown  that  the  limiting  curve  defined  by  Eq,  (197) 
is  the  envelope  of  the  frequency  contours  in  this  region  and  represents 
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the  locus  for  which  the  two  characteristic  frequencies  uj.  and  U)  are 

l  Zr 

equal.  The  frequency  values  along  this  curve  are  indicated  in  Fig.  5 
DISCUSSION 

The  mapping  of  the  R  ,  R.  plane  presented  in  Figs.  2  through  5  en 

x  y 

ables  the  user  to  determine  whether  or  not  a  given  configuration  is 
stable  and  if  stable,  to  determine  the  characteristic  pitch,  yaw  and 
roll  frequencies. 
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IO.  abstract 

An  investigation  of  several  sources  of 
perturbation  on  the  performance  of  a 
gravity  gradient  stabil ‘ nation  system, 
including  mirro-meteoro*  impact,  solar 
radiation  pressure,  station-keeping 
propulsion,  orbital  eccentricity,  and  on¬ 
board  rotating  machinery.  Formulas  are 
developed  which  determine  the  magnitude 
of  pitch,  roll,  and  yaw  perturbations 
from  each  of  these  effects.  These  pertur¬ 
bations  do  -ot  necessarily  rule  out 
gravity  gradient  stabilization,  but  their 
magnitudes  for  a  given  vehicle  configur¬ 
ation  should  be  considered  in  the  light 
of  mission  requirements. 
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